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摘要 : 利用 电流 与 空间 磁场 的 对 应 关系 ， 将 电弧 视 为 具有 变 弧 径 多 条 电流 线 集合 。 
基于 毕 奥 陕 伐 尔 定 律 ， 通 过 正 演 分 析 获 取 组 合 电弧 电流 线 在 空间 的 磁场 分 布 。 为 简化 
问题 复杂 度 ， 考 虑 电弧 未 入 栅 片 的 阶段 ， 利 用 多 种 群 遗传 算法 优化 目标 函数 求解 电磁 
逆 回 题 非 线性 方程 组 ， 以 确定 电弧 空间 位 置 。 在 上 述 基 础 上 ， 为 解决 逆 问 题 不 适 定性 ， 
采用 截断 奇异 值 分 解法 反 演 电弧 电流 分 布 ， 并 通过 GCV 准则 获取 正则 化 参数 。 反 演 结 
果 有 助 于 理解 开关 电弧 的 物理 特性 。 

关键 词 : 开关 电弧 反 演 不适 定性 ”多 种 群 遗 传 算法 ”截断 奇异 值 分 解 ”GCV 准则 
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Truncated Singular Value Decomposition 
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Abstract: By using the corresponding relationship between the current and the space 
magnetic field, the arc 1s regarded as a set of multiple current lines with variable arc 
diameters. Based on Biot-Savart's law, the magnetic field distribution of combined arc 
streamlines in space is obtained by forward analysis. To simplify the complexity of the 
problem, the research only considers the phase of the arc not entering the splitter plates, 
and the multi-population genetic algorithm(MPGA) is used to optimize the objective 
function for solving the nonlinear equations of the electromagnetic inverse problem to 
determine the space position of the arc. On the basis of the above, the truncated singular 
value decomposition(TSVD) method is used to inverse the arc current distribution, with 
the regularization parameters obtained by the GCV criterion. The inversion results will 
contribute to understand the physical characteristics of the switch arc. 
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1 引言 


断路 器 作为 电力 系统 中 的 控制 与 保护 装置 ， 在 
故障 情况 下 切断 线路 以 保护 人 身 和 系统 安全 。 作 为 
有 触 点 开关 ， 断 路 器 在 切断 电流 时 伴随 触 头 分 离 ， 
不 可 避免 地 产生 电弧 中。 电弧 的 燃烧 是 在 多 物理 场 
耦合 作用 下 复杂 的 非 线性 过 程 ， 有 具有 随机 性 和 混沌 
性 品 。 常 用 低压 断路 器 ， 是 以 空气 为 灭 弧 和 绝缘 介 
质 ， 其 电弧 一 方面 给 系统 电路 提供 导电 通路 ， 男 一 
方面 电弧 将 烧 蚀 电极 和 灭 弧 系统 内 结构 部 件 ， 影 响 
开关 寿命 中。 电弧 行为 时 空域 演化 研究 ， 对 开关 电 
器 开 断 能 力 和 寿命 有 着 重要 影响 。 

目前 ， 开 关 电 弧 实 验 研 究 主 要 手段 有 高 速 电荷 
耦合 器 件 (Charge Coupled Device, CCD) 摄像 所、 
光纤 阵列 局、 光谱 测试 号 以 及 磁 测 试 方法 等 。 其 
中 ， 基 于 光学 的 测试 方法 会 影响 实验 样机 的 完整 性 ， 
进而 影响 灭 弧 室 本 体内 电弧 等 离子 体 的 动态 特性 。 
磁 测 试 技 术 多 采用 非 接 触 式 测量 手段 ， 基 于 电弧 电 
流 所 产生 的 磁场 时 空 分 布 反 演 电弧 电流 可 能 存在 的 
空间 位 置 与 形态 "文献 [11-12] 对 比 了 各 优化 算 
法 重 构 电弧 的 误差 ， 采 用 模拟 退火 算法 较 准 确 地 歼 
取 电 弧 形 态 ， 并 利用 磁 成 像 测量 装置 检验 算法 的 可 
行 性 ， 文献 [13-14] 将 电弧 发 生 区 域 痢 分 ， 并 将 电 
磁 正 演算 子 离散 化 为 线性 算 子 ， 采 用 正则 化 方法 反 
演 电弧 电流 ; 文献 [15] 比较 了 各 正则 化 参数 选取 策 
略 重 构 低 压 开关 电弧 电流 的 精度 ， 并 分 析 了 采样 距 
离 以 及 采样 点 数量 对 反 演 精 度 的 影响 。 

综 上 所 述 ， 低 压 开 关 电 弧 模 型 假设 为 变 弧 径 多 
电流 线 集成 ， 基 于 电弧 电流 与 磁场 的 对 应 关系 ， 将 
求解 非 线性 方程 组 转化 为 优化 问题 ， 采 用 多 种 群 遗 
传 算 法 获取 电弧 位 置 及 形态 分 布 以 建立 正 演算 子 。 
为 克服 反 演 的 不 适 定性 ， 采 用 截断 奇异 值 分 解 并 利 
用 GCV 准则 求 取 截 断 参 数 ， 借 助 于 磁场 分 布 捕捉 
电流 分 布 。 结 果 证 实 了 该 方法 能 较 准 确 地 重 构 电弧 
电流 ， 有 助 于 指导 高 性 能 可 靠 性 开关 电器 的 设计 。 


2 ”低压 开关 电弧 物理 模型 


2.1 模型 与 相关 假设 

低压 断路 器 灭 弧 系统 二 维 简 化 模型 如 图 1 所 
示 。 当 动静 触 头 分 离 时 ， 断 口 间 产生 电弧 ， 随 着 行 
程 增加 ， 电 弧 被 拉 长 ， 直 到 完全 打开 ， 动 触 头 到 达 
最 大 开 距 。 在 吹 弧 气流 的 作用 下 ， 电 弧 运 动 至 灭 
弧 室 某 位 置 。 根 据 文献 [7]， 由 于 炎 弧 室 器 壁 的 限 
制 ， 设 电弧 运动 轨迹 在 x0z 平 面 内 。 电 弧 模 型 简化 
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为 多 条 电流 线 集合 ， 各 电流 为 三 段 式 折线 ， 各 条 折 
线 的 节点 位 于 电极 间距 的 等 分 线 上 。 由 于 近 极 区 的 
电流 密度 大 ， 因 此 对 弧 柱 的 收缩 压力 较 强 ， 故 可 暂 
不 考虑 近 极 区 的 弧 柱 半径 。 若 各 节点 坐标 为 4(Cx。， 
0, z). B(x,, 0, z). Cs 0, z). DG, 0, z). 
如 图 2 所 示 ， 则 电极 间距 确定 时 ， 各 节点 的 z 坐标 
为 已 知 量 。 根 据 文献 [ 9]， 磁 场 采 样 平面 距离 电弧 
发 生 面 大 于 15mm 时 ， 栅 片 等 磁性 材料 对 电弧 反 演 
精度 的 误差 小 于 5%， 因 此 采样 距离 为 15mm。 为 
了 验证 该 方法 的 可 行 性 ， 本 文 仅 研究 电弧 未 入 栅 片 
的 阶段 ， 即 图 1 中 虚线 包含 区 域 RQ， 该 段 区 间 面 积 
为 50mm x 30mm。 磁 采样 点 采用 4 x 6 的 阵列 分 布 ， 
如 图 2 所 示 。 


1 一 静 触 点 “2 一 动 触 头 “3 一 灭 弧 室 器 壁 ”4 一 弧 柱 区 5 一 棚 片 
图 1 电弧 电流 线 模型 


Fig.1 The arc current line model 


L3 E * * * * 
X x XBox x X 磁 采 样 点 | 
0mm 
x x Cx x x x 
p | 
ia Lol x Jie Lol 
| 9—— 50mm 一 一 一 "| 


图 2 磁 采 样 点 分 布 


Fig.2 Distribution of the magnetic sampling points 
2.2 电流 线 元 的 磁场 


根据 毕 奥 萨 伐 尔 定律 ， 空 间 电 流 元 产生 的 磁 感 
应 强度 B 为 


Uo. 
B- cos O -cos 0. 1 
And, i 2) (1) 


ÄP, d, 0, 0; 如 图 3 所 示 , uo 为 真空 磁 导 率 。 
图 3 中， 设 空间 电流 线 48 在 x0z 平 面 内 ,5 点 位 
于 距离 AB 为 a 的 平面 BFGH W. 0, 和 6. 分 别 为 端 
A AX 0, Z) Bæn 0, 27) 与 该 处 到 采样 点 S(x,， 
0, z) 的 夹 角 。 推 必得 S 点 处 B.、B,、B. 分 别 为 


Hol b 
也 ,= — (cos 0, -cos 0, 
“gnd, Bp 1 2) (2) 
Hol b, 
B, 2 ———— (cos 0, - cos 0 3 
^ ad, pt | 2) (3) 
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图 3 电弧 电流 元 “4B” 
Fig.3 The arc current unit "4B" 


Ltd b 

* 4nd, b 

X (2) 一 式 (4) "m, d — €, bi, bz, bs, b, 

cos Q4, cos 0; 为 关于 端点 4、B 坐标 的 非 线 性 表达 
AUI, A (2) ~ 式 (4) 可 表示 为 


Ax=b (5) 


(cos 0, -cos 0, ) (4) 


EPRA x= [h Db. Doo 风 的 分 布 ， 首 
先 需 要 确定 各 电流 线 市 点 的 坐标 。 当 电弧 总 电流 值 
1 已 知 时 ， 可 先 假设 各 支线 电流 值 一 定 ， 为 =n, 
n 为 电流 线条 数 。 通 过 智能 优化 方法 搜索 电弧 位 置 
即 各 电流 线 坐 标 ， 算 子 4 即 确定 。 再 通过 解 线性 方 
程 组 反 演 x。 由 于 栅 片 等 铁 磁 材料 对 B, 的 扰动 较 小 ， 
H. B, 相对 于 其 他 磁 感 强度 分 量 的 值 较 大 ， 因 此 在 电 
磁 正 演 分 析 中 ， 采 样 点 测试 采用 分 量 B. 的 分 布 。 


3 ”多 种 群 遗传 算法 优化 


x (5) 可 以 转化 为 下 面 的 优化 问题 ， 即 
min [4x-5| (6) 


XX (6) m, Wh. ÆA irik 
五 搜索 算 子 4 的 优化 过 程 涉及 多 参数 及 复杂 运算 表 
达 式 的 这 一 特殊 性 ， 传 统 的 非 线 性 优化 方法 在 解决 
这 类 问题 时 依赖 对 多 变量 求 极 值 点 的 大 量 计算 。 采 
用 启发 式 的 随机 搜索 方法 能 够 克服 这 一 缺点 ， 在 这 
类 方法 中 ， 遗 传 算 法 (Genetic Algorithm, GA) 是 
一 类 借鉴 生物 界 进化 规律 〈 适 者 生存 、 优 胜 劣 汰 
遗传 机 制 ) 演化 而 来 的 随机 化 搜索 方法 ， 是 一 种 有 具 
有 并 行 机 制 的 全 局 寻 优 方法 。 但 标准 遗传 算法 存在 
着 早熟 收敛 问题 ， 即 某 个 体 的 适 定 度 远 超 其 他 个 体 
时 ， 则 种 群 将 被 该 个 体 控制 ， 群 体 进化 停滞 不 前 。 
为 了 提高 算法 求解 的 质量 ， 在 优化 过 程 中 采用 了 多 
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种 群 遗 传 算法 (Multi-Population Genetic Algorithm, 
MPGA)。MPGA 在 GA 的 基础 上 做 了 以 下 改进 

(1) 多 个 种 群 协同 进化 。 各 个 种 群 取 不 同 的 控 
制 参数 ， 在 GA 中 ， 交 又 算 子 决定 了 全 局 搜索 能 
变异 算 子 决定 局 部 搜索 能 力 ， 交 又 概率 和 变异 概率 
是 常数 。 在 MPGA 中 ,不同 种 群 的 交叉 概率 P. 和 
变异 概率 P, 随机 生成 ， 各 不 相同 ， 即 


[^ 20.74 (0.9—0.7) x rand (N ,1) 


7 
P, 20.0014 (0.05 —0.001) x rand (N ,1) uU 


(2) 各 种 群 通过 移民 算 子 将 进化 过 程 中 各 种 群 
之 间 前 一 种 群 最 优 个 体 代 替 后 一 种 群 中 的 最 差 个 体 ， 
实现 各 种 群 之 间 的 信息 交换 。 

(3) 在 进化 的 每 一 代 ， 利 用 人 工 选择 算 子 选 出 
种 群 的 最 优 个 体 存 入 精华 种 群 。 精 华 种 群 不 进行 选 
择 、 交 又 、 变 异 操作 ， 只 进行 更 新 ， 以 确保 最 优 个 
体 不 被 破坏 和 丢失 。 

MPGA 的 具体 算法 流程 如 图 4 所 示 , 1 为 当前 
代数 ， 初 始 化 种 群 个 数 N= 10， 每 个 种 群 的 个 体 数 
M=40， 最 大 进化 代数 7= 800， 采 用 二 进 制 编码 。 
由 于 电弧 磁场 反 演 问题 的 不 适 定性 ， 可 能 导致 各 电 
流 线 的 位 置 交 又 。 为 此 ， 需 对 算法 变量 施加 约束 。 


目标 函数 /x) 
V use 
170, M, 
P(0,1), EUN 
Y 
个 体 适 应 度 评价 
个 体 适 应 度 个 体 适 应 度 个 体 适 应 度 
Fitness(P(t,1),M)| |Fitness(P(t,2),M) Fitness(P(t,N),M) 
Y Y l 
选择 交叉 变异 | “选择 交叉 变异 … | | 选择 交叉 变异 
TERHET TERET 
Y Y 
| | 下 一 代 种 群 | | | 下 一 代 种 妊 下 一 代 种 妊 
t=t+1,P(t,1) t=t+1,P(t,2) t=t+1,P(t,n) 
mi^ Mn 
Y 
最 优 个 体 存 入 A 
精华 种 群 N 
Y 


4 多 种 群 遗 传 算法 流程 
Fig.4 The fow chart of MPGA 
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通过 引入 罚 函 数 ， 降 低 不 可 行 解 适应 度 以 减 小 被 选 
择 概率 ， 将 约束 问题 转化 为 无 约束 求 极 值 问 题 。 


一 般 约束 优化 问题 为 
min f (x) xER” 
st. h(x)=0 i€ E-(L- 人 (8) 
g(x)20  i€Isí(l-m) 
可 行 域 天 记 作 为 


F-(x€R'|h(x)-0( € Eg, (x) 206€1) (9) 


构造 函数 p(x,0)=f(x)+o p(x)， 其 中 f(x) 为 
MARK. PORTIRA, H. 


PG)- Xh (+ [minio GN" (10) 


其 中 ，a= 1 为 惩罚 因子 。 

经 过 500 代 的 优化 过 程 ， MPGA 与 GA 的 目 
标 函 数 进化 曲线 如 图 5 所 示 。 优 化 各 电流 线 节 点 坐 
标的 相对 误差 曲线 如 图 6 所 示 。 可 以 看 出 ，MPGA 


n — — MPGA 
al ---GA 
1 
s 
Ju 
X. 
5 
Ed 0r A 
"A Ss 
S b S 
n -1- s a 
a9 nis J 
-3 idi c B 
E 
-4 


100 200 300 400 500 600 700 800 
进化 代数 
5 目标 函数 进化 过 程 
Fig.5 Evolution process of the objective function 


20 


—a-MPGA 
18r 一 多 一 GA 


相对 误差 (%) 
S 


| € 
P 


0 2 4 


8 L d 
i 
6+ jm a J 
ra * 
4 ] J 
" A A e 
] bs a \ st] 
\ P x x ge 7^ 
B 1 | K 1 
6 8 10 


12 


电流 线 节点 编号 
图 6 MPGA 和 GA 优化 相对 误差 曲线 
Fig.6 The optimization relative error of MPGA and GA 


ChinaXiv 合 作 期 刊 


能 够 有 效 克 服 GA 的 早熟 收敛 。MPGA 优化 各 变 
量 平 均 相 对 误差 为 2.44%，GA 的 最 大 相对 误差 
为 3.92%。 并 且 MPGA 以 较 快 速度 收敛 到 最 优 值 
0.029 0, GA 收敛 到 最 优 值 0.075 3， 说 明 MPGA 能 
够 提高 算法 的 精度 。 


4 ”截断 奇异 值 分 解 反 演 电 弧 电 流 


通过 MPGA 的 智能 优化 过 程 可 以 确定 式 (5) 
中 的 算 子 4。 由 于 优化 得 到 的 算 子 和 原始 算 子 相 比 
存在 一 定 偏差 ， 且 电磁 逆 问 题 是 病态 问题 ， 因 此 
采用 传统 的 解 线性 方程 组 的 方法 反 演 电流 会 导致 
误差 放大 。 首 先 ， 为 了 反映 问题 的 病态 程度 ， 采 
用 离散 Picard 条 件 判 断 。 对 算 子 4 进行 奇异 值 分 
解 ， 即 


A-Youv; (11) 
i=] 


T RU RA | 四 中 相对 于 奇异 值 0; 趋 于 零 的 
速度 更 快 ， 则 方程 满足 离散 Picard 条 件 ， 可 以 通 
过 正则 化 方法 得 到 最 佳 逼 近 解 。 方 程 4x - 的 离散 
Picard 条 件 判定 如 图 7 所 示 。 可 以 看 出 ， 传 里 叶 系 
数 在 后 面部 分 出 现 波 动 ， 说 明 该 处 奇异 值 对 应 分 量 
受 扰动 影响 较 大 ， 需 利用 正则 化 方法 将 该 段 截断 ， 
以 得 到 最 优 解 。 截 断 参 数 也 称 为 正则 化 参数 ， 通 过 
参数 选择 策略 求 取 。 


—@— |u;b| 
10 是 -二 -Cr 
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图 7 EB Picard 条 件 


Fig.7 Discrete Picard condition 


截断 参数 下 对 于 最 终 的 解 起 着 关键 作用 。 当 上 
选择 过 小 时 ， 去 噪 处 理 过 强 ， 解 过 于 平 请 ， 当 大 选 
择 过 大 时 ， 解 包含 了 噪声 影响 的 不 稳定 分 量 。 在 品 
声 水 平 未 知 的 情况 下 ， 常 用 的 选择 截断 参数 的 方法 
有 工 曲 线 法 和 GCV 准则 。 工 曲线 法 是 通过 绘制 正 
则 解 的 残 差 范 数 与 解 范 数 的 对 数 坐 标 ， 即 
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-log|[A,x - b 
p= iog | 4x-t] m 
-is 

曲率 最 大 的 点 即 对 应 截断 参数 。 但 文献 [19] 
指出 工 曲 线 有 时 过 于 光滑 ， 很 难 找 出 曲率 最 大 的 
一 点 。 

GCV 准则 通过 定义 下 面 的 表达 式 ， 寻 找 其 最 小 
值 对 应 的 大 即 为 截断 参数 UU. B 


G= Jax- (13) 
[trace(I1 - AA; X}? 
AE, A4 为 正则 化 算 子 ，trace RIBERI., XF 
截断 奇异 值 分 解法 ，4i 就 是 将 引起 扰动 的 小 奇异 值 
直接 用 0 FCR, dU ALBUJ XXE AL, ENARAK 
达 式 为 


天 T 
u; b 
— V; 


x,- Ajb- V, (14) 
i-l 


采用 GCV 准则 选取 截断 参数 如 图 8 所 示 。 求 
得 截断 参数 上 =4。 将 大 代入 式 (14) 可 解 得 反 演 电 
流 分 布 如 图 9 所 示 。 


10“ GCV function, minimum at k=4 
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图 8 GCV 准则 选取 截断 参数 


Fig.8 Selection of truncated parameter by GCV criterion 
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9 反 演 电流 分 布 


Fig.9 The inversed current distribution 
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为 了 说 明 算 法 的 精度 ， 对 比 最 小 二 乘法 求解 Ax = 
b 结果 xLs。 最 小 二 乘法 是 一 种 常用 的 求解 超 定 方程 
组 的 方法 ， 甚 解 表示 为 


xis 2 (A A) 47b (15) 


通过 式 (15) npn, 24 k NU A 的 秩 Rank(4) 时 ， 
X (15) 即 为 传统 最 小 二 乘法 的 解 。 但 最 小 二 乘法 
容易 受 问 题 的 不 适 定性 影响 使 得 解 的 稳定 性 变 差 。 
由 图 9 可知， 最 小 二 乘 解 xis 远 远 偏离 预 设 电 流 ， 
而 采用 GCV 准则 选取 正则 化 参数 的 截断 奇异 值 分 
解法 反 演 电流 能 够 反映 预 设 电 流 的 分 布 趋势 。 电 弧 
电流 反 演 相对 误差 见 下 表 。 在 表 中 ， 相 对 误差 6, 为 
表 电弧 电流 反 演 相对 误差 


Tab The relative error of arc current inversion 


电流 线 标号 。 最 小 二 乘法 0%) 。 ”截断 奇异 值 分 解法 06) 
1 8.01 1.14 
2 28.99 2.17 
3 30.85 3.28 
4 39.83 6.56 
5 31.89 6.01 
AA (16) 
Xref 


NP, Xe 为 反 演 电 流 ; xiot 为 预 设 电 流 。 可 见 相 比 
于 最 小 二 乘法 ， 截 断 奇 异 值 分 解 对 于 电磁 逆 问 题 的 
病态 性 具有 较 好 的 抗 噪 能力， 重建 电流 精度 较 高 ， 
证 明了 正则 化 方法 处 理 电 弧 反 演 问题 的 有 效 性 。 

正则 化 方法 反 演 电流 的 误差 主要 是 由 智能 优化 
算法 的 预 处 理 造 成 的 。 由 于 算 子 是 病态 的 ， 即 使 微 
弱 的 偏差 也 会 造成 解 远 远 偏离 真实 值 ， 无 法 通过 正 
则 化 方法 求 近似 解 。 因 此 智能 优化 算法 的 结果 对 方 
程 的 不 适 定 程度 有 较 大 影响 。 


5 结论 


针对 低压 断路 器 电弧 反 演 的 不 适 定性 ， 采 用 多 
电弧 电流 线 表 征 变 弧 径 电弧 ， 利 用 电磁 对 应 关系 反 
演 电弧 ， 以 简化 计算 复杂 度 。 通 过 多 种 群 遗 传 算法 
与 截断 奇异 值 分 解 相 结合 的 手段 对 电弧 位 置 与 电流 
分 布 进行 反 演 ， 利 用 GCV 准则 选择 正则 化 参数 。 
结果 说 明 该 方法 对 重 构 电 绝 电 流 这 一 类 多 变量 不 适 
定 逆 问 题 误差 较 小 。 正 则 化 精度 依赖 于 非 线 性 智能 
优化 结果 的 优 劣 。 为 了 增加 方法 的 可 行 性 ， 降 低 栅 
片 等 非 线性 材料 对 反 演 结果 的 影响 将 是 下 一 步 工 作 
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的 重点 。 
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